Ultracold Atoms in a Tunable Optical Kagome Lattice 
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Geometrically frustrated systems with a large degeneracy of low energy states are of central 
interest in condensed-matter physics [T] [2] . The kagome net - a pattern of corner-sharing triangular 
plaquettes - presents a particularly high degree of frustration, reflected in the non-dispersive orbital 
bands. The ground state of the kagome quantum antiferromagnet, proposed to be a quantum 
spin liquid or valence bond solid [3tilQj, remains uncertain despite decades of work. Solid-state 
kagome magnets [111 \T2 \ suffer from significant magnetic disorder or anisotropy that complicates the 
interpretation of experimental results. Here, we realize the kagome geometry in a two-dimensional 
optical superlattice for ultracold ^^Rb atoms. We employ atom optics to characterize the lattice 
as it is tuned between various geometries, including kagome, one-dimensional stripe, and decorated 
triangular lattices, allowing for a sensitive control of frustration. The lattices implemented in this 
work offer a near-ideal realization of a paradigmatic model of many-body quantum physics. 



Ultracold atoms trapped within optical lattices exhibit 
correlated phases, such as Mott insulators of bosons [13] 
and fermions [14 and effective antiferromagnetic order- 
ing in one-dimensional chains [15j, which are analogous to 
those studied in condensed-matter physics. While most 
experiments have been performed in primitive Bravais 
lattices, with a single lattice site per unit cell, recent 
works have explored non-standard optical lattices with 
a few-site basis such as the two-dimensional honeycomb 
[16] and checkerboard [17] geometries, and double- well 
superlattices [TS] [19] . Such lattices add a low-energy or- 
bital degree of freedom, corresponding to the assignment 
of a particle to each site within the unit cell, that can 
yield non-trivial ordering and dynamics. 

The kagome lattice is obtained by eliminating every 
fourth site from a triangular lattice of spacing a/2, with 
the eliminated sites forming a triangular lattice of spac- 
ing a. The remaining sites generate three connected s- 
orbital bands within a bandwidth on the order of the 
intersite tunneling energy. Intriguingly, the frustration 
of antiferromagnetic interactions in the kagome geome- 
try implies that one of these bands be non-dispersing. 
Such flat bands accentuate the role of interparticle inter- 
actions, leading possibly to crystalline ordering [20 and 
supersolidity [21 for scalar bosons, and ferromagnetism 
of itinerant fermions [22] . 

Here, we realize a two-dimensional kagome lattice for 
ultracold atoms by overlaying two commensurate trian- 
gular optical lattices generated by light at the wave- 
lengths of 532 and 1064 nm (Fig.[l}. Stabili zing and tun- 
ing the relative position of the two lattices, we explore 
different lattice geometries including a kagome lattice, 
a one-dimensional (ID) stripe lattice, and a decorated 
triangular lattice. We characterize these geometries us- 
ing Kapitza-Dirac (KD) diffraction and by analyzing the 
Bloch-state composition of a superfluid released suddenly 
from the lattice. The Bloch-state analysis also allows us 



to determine the ground-state distribution within the su- 
perlattice unit cell as the potential is tuned between ge- 
ometries. The atom-optical tools developed in this work 
are useful for the characterization of other superlattices 
and for superlattice-based atom interferometry. 

Three plane waves of light of equal intensity /, equal 
frequency, and wavevectors (and linear polarizations) ly- 
ing in a plane and intersecting at equal angles, pro- 
duce a triangular lattice of points with zero intensity, 
and a honeycomb lattice of points with maximum in- 
tensity |/ separated by a triangular lattice of intensity 
saddle points with intensity 4/. In our work, the short- 
wavelength (SW) lattice is formed by light that is blue- 
detuned from the principal atomic resonances of rubid- 
ium, so that atoms are attracted to the triangular lat- 
tice of zero-intensity sites with a lattice spacing of a/2 = 
(2/3) X 532 nm = 355 nm. The long- wavelength (LW) lat- 
tice is red detuned, so that its zero-intensity points are 
potential-energy maxima for rubidium atoms. A unit cell 
of the LW lattice contains four sites of the SW triangular 
lattice, labeled A,B,C and D in Fig. [Tj Aligning the po- 
sitions of the LW potential maxima to coincide with sites 
D lowers the potential energies Va,b,c at the indicated 
sites by equal amounts AV = Vd — Va,b,c = § Vlw 
where Vlw is the maximum scalar potential depth of the 
LW lattice (we ignore the ~ 1% vector shift in this lattice 
[T6]). As AV is increased, atoms are excluded from sites 
D, while the remaining sites form the kagome optical lat- 
tice. The kagome geometry persists until Vlw > ^Vsw^ 
at which point atoms become preferentially confined in 
the LW honeycomb lattice. 

Compared with previous proposals [23l [24], our sim- 
pler approach to creating a kagome lattice allows one to 
tune the lattice geometry, thereby controlling its degree 
of frustration. Aligning the LW potential maxima with 
the SW lattice saddle points disfavors population in two 
sites of the four-site unit cell (e.g. Vb^c < ^a,d) pro- 
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FIG. 1: Tunable optical kagome lattice. Three bichromatic light beams intersecting at 120° form a kagome optical lattice 
for ultracold ^^Rb atoms, with the two-dimensional optical potential V^(r) shown in (a). Profiles of the potential of the SW, 
LW, and combined lattices are shown in (b). Sites D of the SW lattice are emptied as their energy offset AV^ exceeds the 
gas chemical potential, so that the remaining sites A, B and C form the kagome geometry, (c) Different lattice geometries 
are created for intermediate LW lattice depths (Vlw < ^Vsw)) by displacing the potential maxima of the SW lattice to the 
high-symmetry points X, Y or Z within the unit cell, while, for higher LW lattice depths, a honeycomb geometry prevails. 



ducing a ID stripe lattice (Fig. [T]^). Aligning the LW 
potential maxima with the SW potential maxima dis- 
favors population in three sites of the four-site unit cell 
(e.g. Va,b,d > Vc)^ producing a decorated triangular lat- 
tice with lowest-energy sites forming a triangular lattice 
while the remaining sites form a kagome lattice of local 
potential minima. 

Experiments were conducted with scalar Bose-Einstein 
condensates (BECs) of ~ 3 x 10^ ^^Rb atoms produced at 
temperatures of about 80 nK in a red-detuned crossed op- 
tical dipole trap with trap frequencies of {ux^oOy^Uz) = 
27r X (60, 30, 350) Hz, with uj^ applying vertically. The 
relative position of the LW and SW lattices was mea- 
sured using two separate two-color interferometers and 
stabilized actively by displacing mirrors in the relevant 
optical paths. The relative position was tuned along two 
linearly independent directions by adding within each in- 
terferometer a variable relative phase shift between the 
two lattice colors (see Methods). 

The first atom-optical method we used to demon- 
strate the tunable superlattice is Kapitza-Dirac diffrac- 
tion [25l [26] , for which the lattice potential is suddenly 
switched on and off, after which the condensate is re- 
leased from the optical trap and imaged after a time of 
flight to evaluate its momentum-space distribution. The 
pulse is sufficiently short that one can neglect the kinetic 
energy of the diffracted momentum components over the 
pulse duration r. The effect of the pulse is then to im- 
print a phase —V{r)r/h proportional to the lattice po- 
tential V{r) onto the condensate wavefunction (which is 
initially nearly uniform). 

The corresponding momentum-space distribution is 
sensitive to the relative displacement of the LW and SW 
lattices. To exhibit this sensitivity we blocked one of 



the incident bi-chromatic lattice beams and examined 
the resulting one-dimensional superlattice formed by two 
bichromatic light beams intersecting at an angle of 120°, 
with potential energy given as V{x) = VLiysin^(g(x + 
Sx)/2) — l/s-vt^ sin^(gx) where 27r/q = 614 nm is the ID 
LW lattice spacing, and Sx is the distance between the 
LW and SW intensity minima. Considering diffraction 
up to second order in the phase modulation depth, the 
atomic populations at wavevectors are given as 

where Jn is the n^^-order Bessel function, evaluated at 
4^LW,sw = yLW.sw ^ /'^^ according to the superscript. 
The lack of inversion symmetry of the lattice produced 
by an incommensurate value of 5x appears as a left/right 
momentum asymmetry in the diffracted matter wave 
(Fig.i. 

A second method to characterize the optical superlat- 
tice is the momentum- space analysis of a superfiuid occu- 
pying the ground state of the lattice potential. For such 
analysis, the optical lattice potential depth was ramped 
up from zero over 90 ms, held at a constant level for 
100 ms, and then suddenly switched off. Imaging after 
a time of flight reveals coherent momentum peaks of the 
condensate wavefunction at the reciprocal lattice vectors 
G of the LW lattice. 

Varying the relative position of the two lattices we 
identify the three high-symmetry lattice configurations 
(Fig. [3|i). Given that the scalar condensate occupies the 
ground state of the lattice potential, its wavefunction can 
be taken as real and positive; thus, its momentum dis- 
tribution is symmetric under inversion. Expansion from 
both the kagome and the decorated triangular lattices 
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FIG. 2: Kapitza-Dirac diffraction from a ID super- 
lattice. Atom diffraction patterns, formed by r — 8 fis 
pulse of the lattice potential followed by 26 ms time of 
flight, exhibit left /right momentum asymmetry (defined as 
(P+g — P_g)/(P+g + P-q) ^^at varlcs as the displacement Sx 
between the LW- and SW-lattice intensity minima is varied, 
reflecting broken inversion symmetry of the lattice, in close 
agreement with the predicted behavior (solid line). 
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shows the three-fold rotational symmetry of the opti- 
cal superlattice. In the ID stripe geometry, one expects 
equally weak occupation of two sites (e.g., A and D), 
and equally strong occupation of the other two sites (B 
and C) of the superlattice unit cell. Such a distribution 
is (nearly) invariant under displacements of a/2 along 
the A-D axis, and condensate diffraction along that axis 
should reflect the shorter periodicity of the SW lattice. 
The momentum distribution should also be symmetric 
under reflection about the A-D axis. Both traits are ob- 
served experimentally. 

The Bloch-state momentum distributions allow one to 
quantify the ground-state wavefunction within a unit 
cell of the superlattice, which we express as ?/;(r) = 
XIq, '0q!'^q;(^ — Sa) where Wa{r) is the normalized Wan- 
nier state wavefunction, the position and IV^^P the 
fractional atomic population of site a G {A, 5, C, D} of 
the unit cell. At low Vlw /Vsw ^ we approximate = w 
as cylindrically symmetric and identical for all a. From 
the momentum-space populations Pg. {i G {1,2,3}) in 
the three first-order diffraction peaks of the LW lattice 
- corresponding to the inner hexagon of peaks in time- 
of- flight images - and that at zero wavevector Pq, one 
determines the distinct quantities 

2Po \w{G,W 

(2) 

where w{0) and w{Gi) are now Fourier components of 



FIG. 3: The momentum-space composition of a BEC 
for various lattices, (a) The kagome and decorated trian- 
gular lattices maintain three-fold rotational symmetry in con- 
figuration and momentum space, while the symmetry of the 
ID stripe lattice is reduced to a parity symmetry (left-right 
in the images). For each setting, a schematic distinguishes 
between sites of high (green) and low (red) atomic popula- 
tion. The lattice potential and expected momentum distri- 
bution are calculated from measured values of the LW and 
SW lattice depths. The area of the white dot reflects the mo- 
mentum population in the expected diffraction pattern, (b) 
Translating the LW-lattice potential maxima along the axis 
connecting sites A and D tunes the lattice between kagome 
and ID stripe geometries, as revealed by the population ra- 
tios Pi identified by color according to the inset. The data 
(averages over 4-5 repeated measurements) agree with calcu- 
lations of the single-particle ground state (solid lines) with 
AV — 14 kHz and Vsw — 40 kHz measured by Kapitza-Dirac 
scattering from LW-only or SW-only lattices. AV was higher 
than the chemical potential /j ^ hx3.5 kHz of the condensate 
in the SW-only lattice, so that modifications of the transverse 
condensate wavefunction due to interactions are negligible. 



the Wannier function, and /3, 7, (5 and e label the four 
sites so that • (s/3 — s^) = 0. Together with the nor- 
malization I'^aP = 1 these quantities determine the 
atomic distribution in the unit cell. We note that invert- 
ing Eqs. |2] gives a discrete degeneracy of solutions that 
can be resolved by prior knowledge of the approximate 
populations. 

We measured the population ratios Pi as the superlat- 
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FIG. 4: Momentum population ratios and atomic dis- 
tributions within the unit cell for a kagome lattice. 

The superlattice was converted from a SW triangular to a 
kagome lattice by increasing Vlw- The D-site population is 
extinguished (b) as the energy offset Ay exceeds the conden- 
sate chemical potential (// ^ /i x 3.5 kHz), and the momentum 
population ratios reach the asymptotic value of 1/9 expected 
for a kagome lattice (a). Data points represent averages of 7 
-10 measurements. Dashed curves indicate expected popu- 
lation ratios Pi for the ideal kagome configuration while the 
shaded region indicates the expected variation in Pi given a 
shot-to-shot instability of ^20 nm in the relative position of 
the LW and SW lattices. 



tice geometry was gradually tuned, using data from the 
second-order diffraction peaks to characterize the Wan- 
nier state (see Methods). Translating the relative po- 
sition of the two lattices (Fig. [s})), one advances from 
the kagome geometry, with equal population in the three 
ratios, to the ID stripe geometry, with two identically 
small ratios, and then to another kagome-geometry lat- 
tice. Our data agree with a calculation of the ground- 
state single-particle wavefunction for the known lattice 
depths. 

We focus finally on the kagome-geometry lattice align- 
ment, and examine the transition between the trian- 
gular and kagome geometries (Fig. |4|. At zero LW 
lattice depth, the atoms are confined in a SW trian- 
gular lattice, and the first-order LW lattice diffraction 
orders are absent, indicating a unit-cell population of 
{A,B,C,D) = (1/4,1/4,1/4,1/4). As the LW lattice 
depth is increased, the population ratios Pi increase and 
the kagome geometry is achieved by gradually expelling 
atoms from one site of the unit cell. The population 
ratios tend toward a limiting value of 1/9 that is a hall- 
mark of diffraction from a kagome lattice wherein the 
atoms are distributed as (A, B^C^D) = (1/3, 1/3, 1/3, 0) 
(see Methods). 



Here, the ground state of the kagome lattice does not 
suffer from frustration, allowing for the stability of a 
Bose-Einstein condensate in the lattice. In the future, 
effects of frustration may be explored by transferring 
bosons into the excited s-orbital fiat band, or by chang- 
ing the sign of the hopping energy [27 so that the flat 
band becomes the lowest-energy band. We note that 
the present choice of wavelengths also yields kagome lat- 
tices for the fermionic isotopes ^Li and ^^K. Introducing 
fermions into the lattice at the appropriate fillings would 
place the Fermi energy within the flat band, allowing for 
studies of flat-band ferromagnetism due to repulsive in- 
teractions [22] or of enhanced Cooper pairing for attrac- 
tive interactions [28]. Furthermore, remarkable versatil- 
ity of the 2D superlattice opens new possibilities to em- 
ulate a distorted kagome lattice as well as a nearly ideal 
one enabling us to explore the possible various quantum 
ground states of the kagome quantum antiferromagnet. 



METHODS 

Optical setup Lattice beams with Gaussian spatial 
profiles, lying and polarized in the horizontal plane, were 
centered onto the condensate. The large ~ 100 /im beam- 
waist diameters of the lattice beams ensured that the 
trapping frequencies were not strongly modified by the 
lattice potential. Laser alignment and relative intensi- 
ties were tuned to produce six-fold symmetric diffraction 
patterns of condensates released from LW- and SW-only 
lattices. In the combined lattice, the relative displace- 
ment of the LW and SW lattices was measured using two 
two-color Mach-Zehnder interferometers, one for beams 
1 and 2 and the other for beams 1 and 3, and stabi- 
lizing all four interference signals by feedback to sepa- 
rate piezo-actuated mirrors in the optical paths. A ro- 
tatable glass plate within each interferometer introduces 
a controlled relative shift between the two lattice colors 
that, following stabilization, is then imparted onto the 
lattice beams at the location of the optical lattice. For 
the Kaptiza-Dirac diffraction experiments, the active sta- 
bilization was engaged with the lattice light at negligible 
intensity before the brief high-intensity pulse of the lat- 
tice potential. For the Bloch-state analysis experiments, 
the stabilization was involved at the beginning of the 
ramp-on of the lattice beams. To reduce deleterious ef- 
fects of atomic interactions in the time-of- flight analysis, 
the vertical optical confinement was switched off about 
0.5 ms before all other optical trapping beams, including 
the lattice beams, were extinguished. 

The optical frequencies of the SW and LW lattice 
beams were stable within a few tens of MHz, at the wave- 
lengths of 532.15 nm and 1063.96 nm respectively. While 
the wavelengths are not perfectly commensurate, we note 
that the relative displacement of the SW and LW lattices 
varies by only (18,9) nm over the {dx^dy) = (28, 14) jam 
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Lattice geometry 


Momentum population ratio Pi 


Triangular (SW) 





kagome 


1/9 


ID stripe 


1 


Decorated triangular 


1 


Honeycomb 


1/4 



TABLE L The momentum population ratios for various lat- 
tice geometries assuming the Wannier wavefunctions at all 
sites to be ^-functions in real space. The hallmark of the 
kagome lattice is the limiting value of 1/9 as observed in Fig. 
[4] Although the ID stripe and the decorated triangular lat- 
tices have the same momentum ratios in this table, they can 
be easily distinguished through their rotational symmetries, 
as can be seen in Fig. [3] 



transverse diameter of the condensates in our experi- 
ments. 

Bloch-state analysis Momentum populations at 
the reciprocal lattice wavevectors are measured by in- 
tegrating under the narrow atomic distributions ob- 
served in time of flight images, thereby accounting for 
effects of finite temperature, atomic interactions, and fi- 
nite size during the expansion of the gas. The Wan- 
nier state Fourier components that appear in Eq. [2] are 
determined, assuming the Wannier state to be Gaus- 
sian, from the second-order diffraction populations as 
|^(0)|V|^(G,)P = V2^o/(P2G, +P-2gJ. The limit- 
ing values of the momentum population ratios enable us 
to distinguish different lattice geometries as summarized 
in Table H. 
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